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Upper and lower bounds are derived for the ground-state energy of neutral atoms which 
for Z — > oo both involve the limits of exact Greens functions with one-body potentials. The 
limits of both bounds are shown to coincide with the Thomas-Fermi ground-state energy. 



SI. Introduction 

(N ■ 

> . 

If^ , A very remarkable property of atoms is that in the hmit Z ^ oo, the Thomas- 

[ Fermi energj^'^^'l^ becomes exact.'^''^® Unfortunately, the very ingenious proofs 

of this beautiful result are somewhat complex. We have strived in developing a 
relatively easier, but rather formal, derivation of this fundamental result for neutral 
I atoms by using, in the process, the Greens function corresponding to the Thomas- 

■ Fermi potential. The derivation rests on the fact that elementary scaling properties 
'. of integrals of the Greens function allow one readily to consider the Z ^ oo limit with 

no difficulty. The basic idea is that integrals of the Greens function for coincident 

■ space points involved in the analysis have particularly simple power law behaviour 
i for large Z. This is spelled out in the text. 

qh' For the Hamiltonian of neutral atoms we choose 

g=f fjj-^yf , , (1) 

I We derive upper and lower bounds on the exact ground-state energy of ([T|), which for 

Z ^ oo are the limits of expressions involving integrals of the exact Greens functions 
with one-body potentials. The limits of both bounds are shown to coincide with the 
ground-state Thomas-Fermi energy, thus establishing the result. 

§2. The upper bound 

We consider first the seemingly unrelated problem of a one-body potential with 
Hamiltonian 

h = ^ + Vir) (2) 
2m 
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where V{r) is the Thomas-Fermi potential 
V(r) = 



r — r' 



2/3 



2m 



R 

ZV3' 



and n(r) = Z^ptf(R) is the Thomas-Fermi density normahzed as 

d'^r n(r) = Z. 

The Greens function corresponding to ([2]) satisfies the equation 

1. _ |iv2 + y(r)l G±(rt,r'0) = (^^(r - r')5(t), 
ar zm J 

where, with appropriate boundary conditions, 
G±(rt,r'0) = T 



^je(Tt)Go(rr,r'0;y), r = 



We write 



Go(rr,r'0;y) 



d^k 

(2^ 



gik.(r~r') 



2m 



■T + U{r,T, k) 



We readily see that U satisfies the equation 



dU 



+ y__k. VC/ + — (V[/)" + i— v2[/ = 0, 
OT m Zm ' 2m 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



with the boundary condition C^|^_q = 0. We are particularly interested in the integral 

y"d3rGo(rr,rO;y), (9) 

where exp [ik • (r — r')] in ([7]) becomes simply replaced by 1. Under a scaling we 
have r = R/Z^/^, V{r) = Z^^^v(W), where w(R) is independent of Z. Accordingly, 
to study the large Z behaviour, we carry out the change of variables r ^ R and 
simultaneously substitute r = T/Z^^^. Also, with the change of variables k — > K, 
k = Z2/3K, the product kV = K^T in 1^ remains invariant. With these new 
variables, ([8]) becomes 



dU 
df 



+ v 



m 



ZV3 



K • VrU + 



viu = o. 



2mZ2/3 



(Vi?C/)' + i 



2mZ2/3 



(10) 



Let lim^_»oo = f^oo- Then pO]) collapses to —dUoo/dt + v = 0, whose solution 
is Uoo = vT. Hence for Z oo, the expression in ([9]) becomes simply scaled by 
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Z^^Z'^ = Z. [On the other hand, if we carry out the unitary scale transformation 
k ^ K, k = Z^/^K, dSD leads to U = vT + 0(^-2/3), and [h'^k'^T/2m + U] ^ 
[h'^K'^T/2mZ^/^ + vT + 0{Z-^/^)]. The latter, under the subsequent change of 
variables K Z^/^K, leads to [h'^K'^T/2m + vT + 0{Z~^/^)], giving the same 
expression for Q as before, with an overall scaling by Z.] 

Accordingly, we have the following limits for large Z, as readily verified upon 
substitution of vT for [/, Z — > cxd : 

/2 f'^ (It f 
d'r — J —— Go(rr, rO; V) ^ Z j d^R pxf(R) = Z, (11a) 



J 27nJ^^T-i£ dr 



2 /d^R 



d^K 

(2^ 



2m 



2 



+ v(R) 



G|,/-^^-|K| 



+ e^ J d^R J d^R' ptf(R) |^^^;| Ptf(RO- (Hb) 

Here, the factor 2 multiplying the r-integrals is to account for spin. The r-integrals 
project out the negative spectrum of h. 

Equation (|lla|) in particular is of fundamental importance. It states that for 
large Z, the Hamiltonian h, allowing for spin, has Z (orthonormal) eigenvectors cor- 
responding to its negative spectrum. Let gi{r,a), . . . ,gzi'r,cr) denote these eigen- 
vectors for large Z. Define the determinantal (anti-symmetric) function 

4>z{riai, vzaz) = det [ga{Yp,ap)\ . (12) 
V Z ! 

Since such an anti-symmetric function does not necessarily coincide with the 
ground-state function of the Hamiltonian H va. ^ m. question, the expectation value 
(^(pz \H\(j)z) with respect to 4>z in (fT2l) can only overestimate the exact ground-state 
energy Ez of H, or at best be equal to it. 

We rewrite the Hamiltonian in ([1]) equivalently as 




e 

a=l 



where is defined in ([2|) with variables Tq,, p^. 
Accordingly, 

Vim Z-'^'^Ez ^ lim 



lim V {g^ \K\g^) + hm Z-'I^Fz, (14) 

Z^oo ^ — ' Z— >oo 
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where 
Fz 



I _ ,1 riz{ra, ra) n(r') 

a I I 

+ y 2^ / I _ ^nz{rcr,ra)nz{r a ,r a ) - \nz[ra,r a )| 



nz{ra,r'a') = ^ ^^(r, cj) 5* (r, a') 



a=l 



or 



Fz^ -e' 



d^r d^r' 



r — r' 



n 



>') \^nz{ra,ra)j 

-\ {^nz{ra,ra)^ (^^^ nz(rV, rV) j 



However, we also have 



dr 



2 

hm Z-2y-nz(rCT,rCT) = hm / Q U.^^Q. y\ 

Z^oo ^ Z^oo 27ri J_ooT — ie 

= ptf(R), 



(15) 



(16) 



(17) 



(18) 



hm y2 {g^ \ha\g^) = hm f 2 ^ 

a=l \ A<0 / 



= hm / d^r — r 

z^oo J 2TnJ_^T-\e 

xi^Go(rr,rO;y), (19) 

where Y1\<q^ ™ '^Yli\<o^ is a sum over all the negative eigenvalues of /i in ([2]), 
allowing for multiplicity but not spin degeneracy. The factor 2 takes the latter into 
account. 

From ([14])-([19]) and (fTlbl) . we finally have 

lim Z-y-Ez ^ ^^^Q^ [ d3R {pM^)f' - t d^R ^ 
z^oo lOvr^m y \r K ,j J ^ 

+ y y d^Rd^R' ptf(R)^^^Ptf(R'), (20) 

and the right-hand side is the coefficient of Z"^/^ of the ground-state Thomas-Fermi 
energy. 
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§3. The lower bound 



Given any arbitrary real and positive function pzi^^), we use the following ele- 
mentary text-book bouncPi : 



z 

E- 

a</3 



Z 
a=l 



,3 Pzjr) 

a r 1 



1 



^ /d^rdVpz(r) 



d'r {pz{r)y 



1/3 



(21) 



Here the real function pz{^) may be chosen to be positive and is otherwise arbitrary 
(i.e., may be chosen at will) to the extent that the integrals on the right-hand side 
of (pT|) exist. We conveniently choose it in such a way that pzi^) Z'^Ptf(R} 
for Z oo, which will then coincide with n(r) used above in Consider the 
Hamiltonian h' = p^/2m -|- V' , where 

V'{r) = -^ + e' 



dV^ 



22 



With pz(r) conveniently chosen, V'{r) may be chosen to be a locally square integrable 
function satisfying V'{r) — > for r — > oo. Let ■i/' be a normalized antisymmetric 
function in (ricji, . . . , rzcrz)- Then (pT|) implies that 



a 

3 



1 Jd'rd\'pz{r)j-^pz{r') 



^l/3^2/3g2 



d^r {pz{r)y 



1/3 



(23) 



Consider the lowest energy E of the Hamiltonian Y2a^'a- '^^e Pauli exclusion 
principle comes to the rescue here.® Concerning the Hamiltonian X^^/Iq, the Z 
"non-interacting" electrons (although each interacts with an external potential V') 
can be put, according to the Pauli exclusion principle, in the lowest energy levels of 
X^a ^'a (allowing for spin degeneracy) if Z is less than the number of such available 
levels. If Z is larger, then the remaining free electrons should have arbitrarily small 
kinetic energies to define the lowest energy of h'^. In either case, E ^ 2 -^^ 
where X]a<o defined as above, is now applied to h' . Accordingly, 

lim Z"^/^ (V-li^lV) > lim Kz, 



Z-*oo 



where 



Z^oa 

oo J Q 

-i Go(rr,rO;l^') 

, T — le OT 



(24) 



^-7/3 £! I d^r dV p^(r)^^p^(r') 
Z J r — r' 



:^l/3^-5/3g2 



d^r (pz(r))^ 



1/3 



(25) 
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and Go(rr, rO;y) is defined as above. Also here we have used the equahty on the 
extreme right-hand side of (jl9p . Since the right-hand side of the inequahty (j24p is 
independent of ^, this inequahty holds with ip corresponding to the ground-state 
function of H as well, i.e., with (ip\H\Tp) corresponding to 



rmn (-i/i 



Ez- 



To the extent that > is arbitrary, we choose it conveniently as 



pz{r) = Z2pTF(R)\/l-e-^"«, 



(26) 



(27) 



where a > is an arbitrary scale parameter. We note that /9tf(R) ~ R^^^^ for 
R ^ 0, and that ptf(R') ~ R^^ for R oo. The factor y^l — exp{—ZaR) ensures 
the integrability of the last integral on the right-hand side of ()25p . We estimate the 
latter for Z — > oo as 



1 



Z2/3 



y"d3Rp2p(R) (l-e-^°^) 



1/3 



/ d'K^Rpl^iR) 

JaR^l/Z ^ 



( 



1 -e" 



Z2 



l>aR>l/Z 



d=^Rp^F(R-) 



+ ^[ d3Rp2p(R) 



1/3 



(28) 



The second integral on the right-hand side is at worst logarithmic in Z. Hence the last 
term on the right-hand side of (|25p vanishes for Z — > oo. Since — (l — e~^°'^) ^ — 1, 
the second term (with the minus sign) on the right-hand side of (j25p is bounded 
below by 



-y I d^Rd^R^ ptf(R) ,^ ^ Ptf(RO- 



Finally, we note that since 1 — e — > 1 for Z 



oo, and with V' = Z^^^v'z and 
Mmz^oov'z = 'f^(R)) the limit of the first expression on the right-hand side of (I25p 
coincides with that in (jllbp for Z — > oo. All told, we see that the lower bound in 
(|24p coincides with the upper bound in (j20p . This completes our demonstration. 

In a future report, we will investigate to what extent this analysis may be ex- 
tended to other interactions. 
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